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ABSTRACT 

We analyze UV spectra for a large sample of 578 Type 1 Active Galactic Nuclei and derive Eddington 
ratios, L/Ledd, from the bolometric luminosities and emission line widths for each object in the sample. 
The sample spans five orders of magnitude in supermassive black hole (SMBH) mass, seven orders 
of magnitude in luminosity, and a redshift range from < z < 5. We include a sample of 26 low- 
redshift Narrow-Line Seyfert Is (NLSls) for comparative analysis. The NLSls have slightly larger than 
average L/Ledd ratios (and smaller SMBH masses) for their luminosities, but those L/L^dd values are 
still substantially below the average for luminous quasars. A large fraction (27%) of the objects overall 
have L/Ledd > 1, which might be explained by non-spherically symmetric accretion. We find no trend 
between L/L^dd and either redshift or SMBH mass. Composite spectra sorted by L/L^dd show an 
unusual emission-line behavior: nearly constant peak heights and decreasing FWHMs with increasing 
L/Ledd- This is in marked contrast to the emission-line behaviors with luminosity, SMBH mass, and 
FWHM(C IV), which clearly show trends analogous to the Baldwin Effect: decreasing line peaks and 
equivalent widths with increasing luminosity, SMBH mass, and FWHM. The origins of the unusual 
behavior with L/L^dd are not understood, but one implication is that metallicity estimates based on 
emission line ratios involving nitrogen show no trend with L/ L^dd in the composite spectra created from 
different ranges in L/Ledd- The NLSl composite, however, shows a slightly high metallicity for its SMBH 
mass and luminosity. Our earlier work suggests that host galaxy mass, correlated with SMBH mass and 
AGN luminosity, is the fundamental parameter affecting BLR metallicities. Some secondary effect, not 
related to L/L^dd, must be enhancing the metallicities in NLSls. 

Subject headings: galaxies: active — quasars: emission lines — galaxies: formation 



1. INTRODUCTION 

The central engines of quasars, and more generally, ac- 
tive galactic nuclei (AGNs) are believed to be powered by 
supermassive black holes (SMBHs). Two of the funda- 
mental properties of AGNs are the SMBH mass and the 
accretion rate of material onto the SMBH. Several indi- 
rect methods have been devised to estimate SMBH masses. 
One set of these methods assumes that the broad emission- 
line region (BLR) is in gravitational equilibrium with the 
central source, so that the SMBH mass can be estimated 
by applying the virial theorem, Msmbh = rv'^/G, to the 
measured line widths (Peterson 1993; Peterson 1997; Wan- 
del et al. 1999; Kaspi et al. 2000; McLure & Dunlop 2001; 
Vestergaard 2002). In reverberation mapping studies, 
-Rblr, the radial distance between the central source and 
the BLR can be estimated from the lag time between con- 
tinuum variations and the emission-line response (Peter- 
son 1993; Peterson 1997; Wandel et al. 1999; Kaspi et al. 
2000). These reverberation mapping studies have demon- 
strated an observed relation of i?BLR cx A iA(5100 A)'^ ''' 
that can be used to estimate Rblr for AGNs over a wide 
range of rcdshifts (Kaspi et al. 2000; McLure & Dunlop 
2001; Vestergaard 2002; Corbett et al. 2003; Warner et 
al. 2003). Netzer (2003) has argued that the slope is not 
known to an accuracy better than about 0.15. 

A second set of methods is based on the tight correla- 
tion between the masses of SMBHs and the velocity dis- 
persions, a, of their host galaxy spheroidal components 



(Ferrarese & Merritt 2000; Gebhardt et al. 2000; Merritt 
& Ferrarese 2001; Tremaine et al. 2002). However, stel- 
lar velocity dispersions are not easy to measure for AGN 
hosts, especially at high redshifts. Because of this, meth- 
ods have been devised using proxies of the velocity disper- 
sion, such as the width of the narrow emission line [O III] 
A5007 (Nelson 2000; Boroson 2003; Shields et al. 2003) or 
the bulge luminosity, L^uige (Magorrian et al. 1998; Laor 
1998; Wandel 1999). Early studies showed a large scatter, 
as much as two orders of magnitude between SMBH mass 
and Lbuige (Ferrarese & Merritt 2000). However, more re- 
cent studies that carefully model the bulge light profiles 
of disk galaxies and thus obtain more accurate values of 
Lbuige show less scatter in A^smbh — ibuigc, similar to that 
in the Msmbh — f relationship (McLure & Dunlop 2002; 
Erwin et al. 2002; Bettoni et al. 2003). Recently, SMBH 
mass has also been shown to correlate strongly with the 
global structure of bulges and ellipticals, such that more 
centrally concentrated bulges have more massive SMBHs. 
This relationship is as strong as the Msmbh — cr relation- 
ship with comparable scatter (Graham et al. 2001; Erwin 
et al. 2002). 

Once the SMBH mass has been estimated, the Ed- 
dington luminosity can be calculated as L^dd = 1-26 x 
lO^^A/sMBH (Mq) ergs s"^ (e.g., Rees 1984; Peterson 
1997). Eddington luminosity is the limit in which the in- 
ward gravitational force acting on the gas exactly balances 
the outward radiation force induced by electron scattering. 
It can be thought of as the maximum possible luminosity 
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for an object of mass Msmbh that is powered by spheri- 
cal accretion (Peterson 1997). The Eddington luminosity 
can be exceeded if accretion is not spherically symmetric 
(see §5, also Osterbrock 1989; Begelman 2002; Collin et 
al. 2002; Wang 2003). AGN luminosities should be di- 
rectly proportional to the accretion rate, L oc Mace, and 
therefore the ratio, L/Ledd oc Macc/M, is an indirect mea- 
sure of the accretion rate relative to the critical Eddington 
value. 

Narrow Line Seyfert Is (NLSls) are a subclass of Seyfert 
Is that exhibit distinct and unusual properties: very nar- 
row broad emission lines {Hfi FWHM < 2000 km s~^) 
with [O III] A5007 / Hf3 ratios of loss than 3 (to exclude 
Seyfert 2s), strong Fe II emission, and unusually strong 
big blue bumps (Osterbrock & Pogge 1985; Kuraszkiewicz 
et al. 2000; Constantin & Shields 2003). NLSls also 
land at one extreme end of the Boroson & Green (1992) 
Principal Component 1 (PCI). It has been suggested 
that PCI is strongly correlated with L/L^dd (Boroson & 
Green 1992; Boroson 2002; Shemmer & Netzer 2002; Con- 
stantin & Shields 2003). Several studies have suggested 
that NLSls have low SMBH masses for their luminosi- 
ties, and thus very high Eddington ratios, near 1 (Mathur 
2000; Kuraszkiewicz et al. 2000; Shemmer & Netzer 2002; 
Shemmer et al. 2003). It has also been suggested that 
NLSls have unusually high metallicitics for their luminosi- 
ties (see §4.4 and Figure 11 below, Mathur 2000; Shem- 
mer & Netzer 2002; Shemmer et al. 2003). Shemmer & 
Netzer (2002) find that NLSls depart from the nominal 
relationship between metallicity and luminosity in AGNs 
(Hamann & Ferland 1999; Dietrich et al. 2003, in prep), 
with some NLSls indicating metallicitics as high as those 
measured in high- luminosity, high-redshift quasars. Be- 
cause of their high metallicitics and high Eddington ratios, 
Mathur (2000) proposed that NLSls are analogs of high- 
redshift {z > 4) quasars, in that they may both be in an 
early evolutionary phase, residing in young host galaxies. 

We have collected a large sample of 578 spectra of 
"Type 1" AGNs (quasars and Seyfert galaxies with broad 
emission lines) that span the rest-frame UV wavelengths 
needed for this study (Dietrich et al. 2002). We com- 
pute composite spectra from different ranges in the Ed- 
dington ratio, L/Ledd- We include a composite spectrum 
produced from a subsamplc of 26 NLSls for comparative 
analysis. We present measurements of the emission lines in 
these composite spectra and investigate their relationship 

to L/Ledd- 

2. EDDINGTON RATIO DETERMINATIONS 

We first estimate SMBH masses by applying the virial 
theorem, Msmbh = rv^/G, to the line-emitting gas (for 
more details, see also Kaspi et al. 2000; Peterson & Wan- 
del 2000; McLure & Dunlop 2001; Vestergaard 2002; Cor- 
bett et al. 2003; Warner et al. 2003). Kaspi et al. (2000) 
express the SMBH mass as 
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-Rblr is the radial distance between the BLR and the cen- 
tral source, and FWHM applies to the broad emission line 
profile. We estimate i?BLR based on the observed relation 
between Rblr for a particular line and the continuum lu- 
minosity (Wandel et al. 1999; Kaspi et al. 2000; Vester- 
gaard 2002; Corbett et al. 2003). A particular line must 
be specified because reverberation studies have shown that 
the BLR is radially stratified, such that higher ionization 
lines tend to form closer to the central engine than lower 
ionization lines (Peterson 1993). 

We select the C IV A 1549 emission line instead of H(3 to 
estimate SMBH masses because it is more readily observed 
across the entire redshift range from z ^ to z 5. We 
find that there is approximately a 1:1 correlation between 
the SMBH mass obtained from C IV and that obtained 
from Hp. There can be significant deviations from this 
for individual objects, but the relation holds well for av- 
erages of many objects and for measurements based on 
composite spectra (see Vestergaard 2002; Warner et al. 
2003 for further discussion). 

Reverberation studies indicate that the radius of the 
BLR for C IV is about half that of Hp (Stirpe et al. 1994; 
Korista et al. 1995; Peterson 1997; Peterson & Wandel 
1999). We therefore modify the equation given by Kaspi 
et al. (2000) to obtain 
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Sec Warner et al. (2003) for more details. From Equations 
(1) and (2), we derive 
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set of cosmological parameters and those suggested by WMAP results 
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Vestergaard (2002) calibrated mass derivations based on 

C IV against estimates using FWHM(iJ/3) and direct mea- 
surements of Rblr{H(3) from reverberation mapping. Her 
technique yields essentially the same mass relationship 
(within 10%) as Equation 3, which helps to confirm the 
factor of 2 scaling adopted here between i?BLR(C IV) and 
Rblr{HP). Vestergaard (2002) finds that SMBH masses 
estimated by applying this equation to single-epoch spec- 
tra of individual objects have a Icr uncertainty of a factor of 
three when compared to studies that use Hp and a direct, 
reverberation measure of the BLR radius. Our compos- 
ite spectra average over variabilities and object-to-object 
scatter, which should significantly reduce the uncertain- 
ties. See also Krolik (2001), Netzer (2003), Corbett et al. 
(2003), and Vestergaard (2004) for further discussion of 
the uncertainties. 

We next estimate bolometric luminosities, L, based on 
an integration over a typical quasar continuum shape. We 
use the cosmological parameters Hq = 65 km s~^ Mpc~^, 
= 0.3, and Ha = (Carroll, Press, & Turner 1992) 
throughout this paper^. We assume a segmented power- 
law of the form oc j^" to approximate the continuum 
shape, with a = —0.9 from 0.1 A to 10 A, q = —1.6 from 
10 A to 1000 A, and a = -0.4 from 1000 A to 100,000 
A (Zheng et al. 1997; Laor et al. 1997; Brotherton et al. 
2001; Vanden Berk et al. 2001; Dietrich et al. 2002). It 

studies in the pre- WMAP era. Furthermore, the difference between our 
in a difference of less than 20% in luminosity for a wide redshift range 
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is now well established that the mean UV-IR slope is a 
function of luminosity (see Dietrich et al. 2002). For the 
luminosity range spanned by our sample, the average pow- 
erlaw index between 1000 A and 100,000 A ranges from ~ 
-0.2 to ~ -0.6 (M. Dietrich, private communication). This 
range could cause a scatter of ~ 20% in our estimates of 
L. 

Integrating this nominal spectrum over all wavelengths 
implies bolometric corrections of 4.36 and 9.27 for 
ALa(1450A) and ALa(5100A), respectively. These cor- 
rections are slightly lower than the correction of 11.8 to 
ALa(5100A) derived by Elvis et al. (1994), but in good 
agreement with more recent derivations (Kaspi et al. 2000; 
Vestergaard 2004). We use this bolometric correction and 
SMBH masses from Equation (3) to obtain Eddington ra- 
tios: 
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3. DATA & ANALYSIS 



Our sample is comprised of 578 Type 1 (broad-line) 
AGN spectra with rest-frame UV wavelength coverage 
that encompasses the range 950 < A < 2050 A. The spec- 
tra were obtained by several groups using various ground- 
based instruments as well as the Hubble Space Telescope 
(HST) and the International Ultraviolet Explorer (lUE) 
(see Dietrich et al. 2002 and 2004, in prep for more de- 
tails). The sample spans a redshift range from ^ 2 < 5, 
seven orders of magnitude in luminosity, and five orders 
of magnitude in SMBH mass. One unique aspect of this 
sample is that it contains new observations of faint quasars 
at redshift z > 2.5 (e.g., Steidel et al. 2002, Dietrich et 
al. 2002). Thus we can avoid to some degree the bias to- 
ward higher luminosities at higher redshifts inherent in 
magnitude-limited samples. The sample spans at least 
three orders of magnitude in luminosity at all redshifts 
(see Fig. 1 in Dietrich et al. 2002). We determined the 
radio loudness for the quasars using the radio flux den- 
sities given in Vcron-Cetty & Veron (2001). We used 
the definition of radio loudness given by Kellermann et 
al. (1989). Our classifications of radio-loud quasars are 
consistent with classifications available in the literature 
(Wills et al. 1995; Bischof & Becker 1997; Wilkes et al. 
1999; Stern et al. 2000). 

We use an automated program to estimate the FWHM 
of C IV in each spectrum (see Warner et al. 2003 for de- 
tails). Comparisons between the FWHMs estimated by 
the program and those measured manually indicate an er- 
ror of < 10% in the automated results. Lines containing 
significant absorption are flagged by the program and their 
FWHMs are estimated manually (by interpolating across 
the absorption feature). We use the FWHM of C IV and 
the continuum luminosity, XL\{1450A) to estimate the 
central SMBH mass and L/L^dd for each quasar based on 
the equations given in §2. 

We then sort the quasars by L/L^dd into seven bins 
(see Figure 1): L/L^ < 0.25, 0.25 < L/Lm < 0.33, 
0.33 < L/L^dd < 0.50, 0.50 < L/L^dd < 0.67, 0.67 < 
L/Ledd < 1.00, 1.00 < L/Ledd < 2.00, and L/L^m > 2.00, 
and compute seven composite spectra. Each composite 
spectrum is the average of all the quasar spectra in a bin. 
Table 1 lists various parameters for the composites, includ- 



ing the mean values of Msmbh, L, FWHM(C IV), L/Ledd, 
and the redshift, z, as measured from the individual ob- 
jects contributing to each composite. Also listed are the 
numbers of objects contributing at the wavelength of the 
C IV emission line. The spectral slopes, a, are measured 
from each composite spectrum and constrained by the flux 
in 20 A wide windows centered at 1450 A and 1990 A. 

Calculating composite spectra significantly improves the 
signal-to-noise ratio and averages over object-to-object 
variations. Since narrow absorption features may affect 
the emission line profiles in composite spectra, we devel- 
oped a method to detect strong narrow absorption fea- 
tures. The contaminated spectral region of the individual 
spectrum is then excluded from the calculation of the com- 
posite spectrum. For more details about creating compos- 
ite spectra, see Brotherton et al. (2001), Vanden Berk et 
al. (2001), Dietrich et al. (2002), and Warner et al. (2003) 

For comparison with the L / L^dd composites, we also cre- 
ate composite spectra for different ranges in SMBH mass 
(10^' - 10^ Mq, 10^ - 10^ Mq, 10* - 10^ Mq, 10^ - 10^° 
Mq, and > 10^° M©; see also Warner et al. 2003), L 
(10^^ - 10^5 ergs s-\ 10^^ - 10^^ ergs s'S 10^^ - 10^^ 
ergs s-\ 10"^ - 10'**^ ergs s^^, and 10''^ - 10*^ ergs s"!), 
and FWHM(C IV) (< 2000 km s'^, 2000 - 4000 km s'^, 
4000 - 6000 km s'^, 6000 - 8000 km s-\ and > 8000 km 
s~^). See Table 1 for additional information. 

We also create a composite spectrum of 26 NLSls that 
were classified by others according to the criteria described 
in §1 (Kuraszkiewicz et al. 2000; Wang & Lu 2001; Con- 
stantin & Shields 2003). This subsample is drawn from 
the same overall sample used to create the other compos- 
ites described above. Eighteen of these objects have data 
at C IV and they span a range in L/Ledd from ~ 0.1 to 2. 

We correct each composite spectrum for strong iron 
emission lines using the empirical Fe emission template 
that was extracted from I Zw 1 by Vestergaard & Wilkes 
(2001), which they very kindly provided for this study (see 
Dietrich et al. 2002 and Warner et al. 2003 for more de- 
tails). The Fe II contribution is generally small at wave- 
lengths < 2000 A, but the correction for this emission 
improves the measurements of weak lines such as N III] 
A1750 and He II A1640. We fit the continuum of each Fe- 
subtracted spectrum with a powerlaw of the form F^, ocv". 
Figure 2 shows the final Fe Il-subtracted composite spec- 
tra normalized by the continuum fits. 

To measure the broad emission lines, we use a spectral 
fitting routine developed in the IDL language, that em- 
ploys minimization. We fit each line with one or more 
Gaussian profiles, with the goal of simply measuring the 
total line strengths free of blends. When necessary, we use 
the profile of strong unblended lines, such as C IV, to con- 
strain the fits to weaker or more blended lines (see Warner 
et al. 2003 for details of our fitting procedure). Figure 3 
shows an example of our fits. 

The continuum location is the primary uncertainty in 
our flux measurements. We estimate the la standard de- 
viation of our measurements of the fluxes of Lya and C IV 
to be <10% based on repeated estimates with the con- 
tinuum drawn at different levels. We estimate the uncer- 
tainty in weaker lines by the same method to be ~ 10-20%. 
There are also secondary uncertainties due to line blend- 
ing, which can be important for some of the weak lines 
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and for N V in the wing of hya. 

4. RESULTS & COMPARISONS 

4.1. L/Lgdd dnd Super-Eddington Accretion 

Figure 4 shows the distribution of L/Ledd for the entire 
sample. A large fraction (27%) of the objects in our sam- 
ple have L/Ledd > 1 (see also Fig. 1). This result is not 
precise because our mass estimates for individual objects 
have factor of ~3 uncertainties (§2), which are comparable 
to the width of the distribution in Figure 4. Nonetheless, 
it is interesting that the sample mean is close to the Ed- 
dington limit, with < L/L^dd > ~ 0.9 (Fig. 4). Also note 
that trends with L/L^dd that we discuss below are more 
reliable than the individual measurements because they 
rely on relative L/L^dd- 

Consistent with previous studies (§1), we find that the 
NLSls have generally high Eddington ratios for their lumi- 
nosities, including several objects with LjL^dd ^ 1- How- 
ever, the NLSls do not have the highest Eddington ratios 
in our sample. Quasars with high luminosities and narrow 
C IV emission lines often have L/L^dd > 2. The most ex- 
treme of these objects, such as BR2248-1242 (see Warner 
et al. 2002), can have derived Eddington ratios approach- 
ing 10 (see Figure 1). 

4.2. Correlations with L/L^dd 

Figure 5 shows the distributions in redshift, bolometric 
luminosity FWHM(C IV), and SMBH mass as a func- 
tion of L/Ledd for the entire sample. L/Ledd correlates 
positively with L and negatively with FWHM(C IV), but 
these correlations may be attributed largely to our deriva- 
tion of L/Ledd from these quantities. In fact, the slopes 
in these correlations are matched well by the parameter 
relationships in Equation 4. In agreement with Woo & 
Urry (2002), we find no trend between L/Ledd and either 
redshift or SMBH mass. We find that the weak trend 
in Figure 5 between L/Ledd and redshift is due to i) a 
trend for larger L/Ledd with increasing L, and ii) a bias 
for more high L objects at higher redshifts in our sample. 
Sub-samples spanning narrow ranges in luminosity show 
that there is no trend between L/Ledd and redshift once 
these biases are removed (see Figure 6). In all four panels 
of Figure 5, there are no clear differences between radio 
loud and radio quiet objects. 

Table 1 shows that there is no apparent trend between 
Eddington ratio and the slope of the UV continuum. The 
NLSl composite, though, exhibits a steeper (softer) UV 
spectrum than the L/Ledd composites (see Table 1). This 
is consistent with findings that NLSls in general have red- 
der spectra than typical Type 1 AGNs (e.g., Crenshaw et 
al. 2002; Constantin & Shields 2003). 

Table 2 lists for each L/Ledd and NLSl composite spec- 
trum the line fiuxes relative to Lya, the rest-frame equiv- 
alent widths (REWs) as measured above the fitted con- 
tinuum, and the FWHMs. Figure 7 plots the REWs of 
selected emission lines as a function of L/Ledd- The NLSl 
composite spectrum is displayed (plotted as a triangle) 
for comparative purposes. The dotted lines are linear fits 
to the L/Ledd composite data (excluding the NLSl com- 
posite). Interestingly, while most emission lines decrease 
in REW with increasing L/Ledd, Lya and O VI exhibit 
a positive trend between REW and L/Ledd, and N III] 



shows no trend at all. O VI has larger measurement errors 
than most other emission lines, so it is unclear whether 
this positive trend between O VI REW and L/Ledd is real 
or not. 

The REWs of emission lines in the NLSl composite (rep- 
resented by a triangle) generally do not match the trend 
fit to the L/Ledd composites in Figure 7. The NLSl com- 
posite is above the fitted trend for some emission lines and 
below the fitted trend for others, regardless of the slopes 
of the trends. 

4.3. Surprising Emission-Line Behaviors 

The composite spectra sorted by L /Ledd show a surpris- 
ing emission-line behavior: nearly constant peak heights 
and decreasing FWHMs with increasing L/Ledd (see Fig- 
ure 2). This is in marked contrast to the emission- 
line behaviors in composite spectra sorted by luminos- 
ity SMBH mass, and FWHM(C IV), which clearly show 
trends analagous to the Baldwin Effect: decreasing line 
peaks and equivalent widths with increasing luminosity, 
SMBH mass, and FWHM. Figure 8 compares the emission- 
line behaviors in these different composites (see also Wills 
et al. 1993; Croom et al. 2002; Dietrich et al. 2002; 
Warner et al. 2003). In particular, the composite spectra 
created from different ranges in FWHM(C IV) clearly show 
a trend analagous to the Baldwin Effect despite spanning 
a range of less than half an order of magnitude in average 
luminosity (Table 1). In contrast, the L/Ledd composites 
span a wide range in FWHM(C IV) and an order of mag- 
nitude in average luminosity but do not show any behavior 
similar to the Baldwin Effect. This suggests that the Bald- 
win Effect may actually be related to SMBH mass (which 
correlates positively with both L and FWHM), since the 
L/Ledd composites have nearly constant Msmbh- 

To illustrate this point further. Figure 9 compares 
L/Ledd composite spectra created for a narrow range in 
SMBH mass (IO^Mq < A^smbh < IO^Mq) and a narrow 
luminosity range (10-*^ ergs/s < L < 10^^ ergs/s). Both 
sets of spectra are shown prior to normalization to the 
continuum. The emission line behavior described above 
(e.g. constant peak heights, etc.) is clearly evident in the 
L/Ledd composites at nearly constant A'/smbh, but not in 
the composites with nearly constant L. The composites at 
nearly constant L show a trend for decreasing peak heights 
and equivalent widths with increasing AfsMBH- 

In composite spectra sorted by luminosity, SMBH mass, 
and FWHM(C IV), the Baldwin Effect is not seen in N V 
(see Figure 8, also Dietrich et al. 2002; Warner et al. 
2003), leading to a higher N V / C IV ratio in objects with 
higher luminosities, SMBH masses, and FWHMs(C IV). 
However, in the L/Ledd composites created from our en- 
tire sample, N V clearly decreases in REW as L/Ledd in- 
creases, yielding a nearly constant N V / C IV ratio across 
the full range of L/Ledd (see Figure 10). Furthermore, 
the L/Ledd composites created from a narrow range in 
AfsMBH exhibit this behavior of nearly constant N V / 
C IV, while the composites created from a narrow range 
in L do not (see Fig. 9). The composites at nearly con- 
stant L show a trend for increasing N V / C IV toward 
lower L/Ledd (higher SMBH masses). The differences in 
emission-line behaviors at nearly constant Msmbh and at 
nearly constant L have implications for the origin of the 
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Baldwin Effect (see Warner, Hamann, & Dietrich 2004, in 
prep). Throughout the rest of this manuscript (including 
figures and tables), "L/Ledd composites" refer to the com- 
posite spectra created from our entire sample and sorted 
by L/Ledd- 

4.4. Metallicities 

We compare emission line flux ratios to plots of metallic- 
ity vs. line ratio based on theoretical models (sec Figure 5 
in Hamann et al. 2002 and Figure 3 in Warner et al. 2002). 
Ratios involving nitrogen lines are especially valuable in es- 
timating metallicity, Z, due to the expected "secondary" 
N production via the CNO cycle of nucleosynthesis in stars 
(Shields 1976; Hamann & Ferland 1992, 1993, 1999; Fer- 
land et al. 1996; Hamann et al. 2002). In the CNO cycle, 
N is produced from existing carbon and oxygen and thus 
the nitrogen abundance scales as N/H oc Z'^ and N/0 oc 
0/H oc Z (Tinsley 1980), providing a sensitive metallicity 
diagnostic even when direct measurements of Z are not 
available (for more discussion, see Wheeler et al. 1989; 
Hamann & Ferland 1999; Henry et al. 2000; Hamann et 
al. 2003; Pilyugin 2003; Pilyugin et al. 2003). We pre- 
fer to base our metallicity estimates on the calculations in 
Hamann et al. (2002) that use a segmented powerlaw for 
the photoionizing continuum shape because this shape is 
a good approximation to the average observed continuum 
in quasars (Zheng et al. 1997; Laor ct al. 1997) and it 
yields intermediate results for line ratios that are sensitive 
to the continuum shape, such as N V/He II. 

Figure 10 shows metallicities inferred from several line 
ratios as a function of L/L^dd- The uncertainties shown in 
Figure 10 derive solely from the la measurement uncer- 
tainties discussed in §3 and do not include the theoretical 
uncertainties in the technique we use to derive metallicities 
from the line ratios. Our best estimate of the overall metal- 
licity from each spectrum (labeled as Average in Figure 10) 
is obtained by averaging the results of the line ratios that 
we believe are most accurately measured and most reliable 
from a theoretical viewpoint. Specifically, we average the 
metallicities derived from N III]/C III], N III]/0 III], and 
N V/C IV (or when available, N V/(C W+O VI)). See 
Hamann et al. (2002) and Warner et al. (2003) for further 
discussion. 

All of the line ratios involving N III] and N V show N/0 

and N/C ratios that are solar or greater. This implies a 
metallicity of > 1 Zq if N is mostly secondary. One impli- 
cation of the unusual emission-line behavior with L/L^dd 
(discussed in §4.3) is that the derived metallicities on av- 
erage show no trend with L/L^dd- The only line ratio to 
exhibit a strong trend with L/L^dd is N V/0 VI. 

Figure 11 shows our best estimates of the overall metal- 
licities (derived from several line ratios, as described above 
for the "Average" in Figure 10) from the composite spec- 
tra sorted by L/L,.dd, SMBH mass, and luminosity. The 
NLSl composite spectrum is plotted in all three panels 
for comparison. The NLSls exhibit a metallicity that is 
slightly high (a 30%-40% enhancement) for their SMBH 
masses and luminosities (see also Shemmcr & Nctzcr 2002; 
Shemmer et al. 2003), but still well below the metallicities 
derived for most luminous quasars. 

It is difficult to compare our results directly to Shemmer 
& Netzer (2002) because they measure only N V/C IV and 



N V/He II ratios for a sample of individual objects. How- 
ever, it does appear that the metallicity enhancement we 
derive as an average from several line ratios is smaller than 
the Shemmer & Netzer result based only on N V/C IV and 
N V/He II. 

5. DISCUSSION 

There is evidence supporting super-Eddington accretion 

rates for ^27% of the objects shown in Figures 1 and 4. 
Previous studies have also found quasars that appear to be 
accreting at super-Eddington rates (e.g., Collin et al. 2002, 
Vestergaard 2004). It is also worth noting that changing 
the index on the Rblr — L relationship (see §1) from 0.7 
to 0.5 (e.g. Nctzcr & Laor 1993; Shields et al. 2003) would 
actually increase our estimates of L/Ledd, by as much as 
a factor of ~ 4 in the brightest quasars (see also Netzer 
2003). Super-Eddington accretion can be explained simply 
by non-spherically symmetric accretion (Osterbrock 1989). 
It has been suggested that accretion disks with radiation- 
driven inhomogeneities could produce luminosities up to 
100 times Eddington (Begelman 2002; Wang 2003). 

Woo & Urry (2002) find that only 9% of their objects (21 
out of 234, see their Fig. 7) have L/Ledd > 1- However, 
this may be due to the selection effect that their sample 
contains only 9 objects with L > 10^^ ergs s~^. In our 
sample, only 13% of the objects with L < 10^^ ergs s~^ 
have L/Ledd > 1 (compared to 8%, 18 out of 225, in Woo 
& Urry 2002), but 41% of the objects in our sample with 
L > 10^'' ergs s~^ have super-Eddington ratios. 

The NLSl composite spectrum exhibits many properties 
unlike the L/Ledd composite spectra. The NLSl compos- 
ite has a much steeper ( "softer" ) UV continuum than the 
L/Ledd composites, consistent with studies finding that 
NLSls generally have redder UV continua than typical 
AGNs (e.g., Crenshaw et al. 2002; Constantin & Shields 
2003). There is no apparent trend between continuum 
shape and L/Ledd (Table 1), but there are trends be- 
tween continuum shape and both SMBH mass (see Table 
1; Warner et al. 2003) and luminosity (Table 1; M. Diet- 
rich, private comm.), such that objects with lower SMBH 
masses and lower luminosities have steeper ("softer") UV 
continua. Therefore, the steeper UV continuum in the 
NLSl composite may be due to the NLSl composite hav- 
ing a much lower average Msmbh and L than the L/Ledd 
composites. 

Interestingly, Figure 2 shows that the NLSl compos- 
ite spectrum most strongly resembles the composite spec- 
trum created from objects with L/L^dd > 2. It seems 
to fit at the top of Figure 2 and not where it would be 
placed based on its average L/Ledd of 0.67. Based on this 
similarity between NLSls and quasars with high L/Ledd, 
it has been suggested that high-redshift (z > 4), high- 
luminosity, narrow-lined quasars are analogs of NLSls 
(Mathur 2000). However, despite the similarity in Fig- 
ure 2 between the NLSl composite and the L/Ledd > 2 
composite, the REWs of some emission lines in the two 
spectra can be quite different (see Figure 7, Constantin & 
Shields 2002). Moreover, the NLSls do not fit the gen- 
eral trends in Figure 7 between emission line REW and 

L/Ledd- 

If host galaxy mass, which correlates strongly with 
SMBH mass (see §1), is the fundamental parameter af- 
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fccting BLR metallicity, it is reasonable to expect no trend 
between metallicity and L/Lf-dd because there is no trend 
between LjL^dd and SMBH mass in the LjL^dd compos- 
ites. The results in Figures 10 and 11 confirm this expecta- 
tion. The different L/ L^dd composites have similar average 
SMBH masses (Table 1) and similar metallicities. How- 
ever, the NLSl behavior is surprising. They have roughly 
the same metallicity as the L/Ledd composites even though 
their average SMBH mass and luminosity are almost two 
orders of magnitude lower. Figure 11 shows more directly 
that the NLSls have slightly high metallicities for their 
luminosities and SMBH masses (see also Shemmer & Net- 
zer 2002; Shemmer et al. 2003). If host galaxy mass, 
correlated with Msmbh, is the main driver behind AGN 
metallicities (Figure 11 and Warner et al. 2003), then 
clearly some other factor is enhancing the metallicities in 
NLSls. The magnitude of the NLSl enhancement is mod- 
est, roughly 30%-40%. Note, in particular, that the NLSl 
metallicities are still well below the values derived for most 
luminous (e.g., high-redshift) quasars. 

Shemmer & Netzer (2002) suggest that high L/Ledd is 
driving the high metallicities in NLSls. However, we have 
shown that there is no correlation between metallicity and 
L/Ledd (Figures 10 and 11). In addition, luminous narrow- 
lined quasars with the highest values of L/L^dd (such as 
BR2248-1242, Warner et al. 2002) do not have high metal- 
licities for their SMBH masses. We conclude that the addi- 
tional parameter affecting NLSl metallicities is not related 

to L/Ledd- 

Constantin & Shields (2003) suggest that the result 
for NLSls having lower metallicities than high-redshift 
quasars runs counter to the hypothesis that these two 
groups of objects are in a similar early evolutionary phase. 
However, because the high-redshift quasars have higher 
SMBH masses and luminosities than NLSls, they would 
naturally be expected to have higher metallicities. If both 
types of objects are in similar early evolutionary phases, 
and/or reside in young or rejuvenated host galaxies and 
the metallicities of both types of objects are thus enhanced 
by a similar amount, then the high-redshift quasars would 
still exhibit a higher metallicity than the NLSls due to 
their higher SMBH masses. Thus, we conclude that more 
information is necessary to determine if either NLSls or 
the high-redshift quasars (or both) are preferentially young 
objects. 

6. SUMMARY & CONCLUSIONS 

We have examined a large sample of 578 AGNs that 

spans five orders of magnitude in SMBH mass, seven or- 
ders of magnitude in luminosity, and a rcdshift range from 
Q < z <f). We estimate SMBH masses using the virial the- 
orem and formulae given in Kaspi et al. (2000), and then 
derive Eddington ratios. To improve the signal-to-noise 
ratio and average over objcct-to-objcct variations, we cal- 
culate composite spectra for different ranges in L/L^dd- 



We include a composite spectrum of a sample of 26 NLSls 
for comparative analysis. Our main results are as follows. 

1) We find that a large fraction (27%) of the objects in 
our sample have L/L^dd > 1- These super-Eddington ra- 
tios may be explained by non-spherically symmetric accre- 
tion. While NLSls generally show high Eddington ratios 
for their luminosities, the objects with the highest L/Ledd 
are high-luminosity, narrow-lined quasars. 

2) There is no trend between L/L^Ad and either red- 
shift or SMBH mass. L/L,.,id does correlate positively 
with luminosity and negatively with FWHM(C IV), but 
these trends may be attributed largely to our derivation 
of L/Ledd from these quantities (see Equation 4). 

3) There is no trend between the shape of the UV con- 
tinuum and L/Ledd- The NLSl composite has a much 
steeper (softer) continuum than the L/Ledd composites. 
This is consistent with a trend between continuum shape 
and SMBH mass. 

4) The composite spectra sorted by L/Ledd exhibit 
an unusual emission-line behavior: nearly constant peak 
heights and decreasing FWHMs with increasing L/Ledd 
(Figure 2). The origins of this behavior are not under- 
stood, but it is in marked contrast to the emission-line be- 
haviors in composite spectra sorted by luminosity, SMBH 
mass, and FWHM(C IV) (Figure 8), all of which clearly 
show trends in the line REWs analagous to the Baldwin 
Effect. 

5) The composite spectra show no trend between L/Ledd 
and metallicity (Figure 10). This is consistent with SMBH 
mass being related to the fundamental parameter affecting 
BLR metalhcity (Warner et al. 2003). 

6) The NLSl composite exhibits several unusual behav- 
iors. It generally docs not fit the trends between emission 
line REWs and L/Ledd as defined by the L/Ledd compos- 
ites. It also has a metallicity that is slightly high for its 
average SMBH mass and luminosity, although still well 
below the high metallicities exhibited by the most lumi- 
nous quasars with the most massive central SMBHs. The 
quasars with the highest L/Ledd^ high-luminosity quasars 
with narrow C IV emission, do not have high metallicities 
for their SMBH masses and luminosities. Our earlier work 
(Warner et al. 2003) is consistent with the theory that 
host galaxy mass, correlated with SMBH mass (and AGN 
luminosity), is the fundamental parameter affecting BLR 
metallicity. We conclude that i) there must be some sec- 
ondary effect cnhandng the metallicity in NLSls, and ii) 
this secondary effect is not related to L/Ledd- 
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Fig. 1. — Distribution of the two measured quantities. The dotted lines represent constant L/Lg^j^ ratios, with values indicated at the 
right. We created composite spectra for each of the seven ranges in Eddington ratio shown above. The filled diamonds represent radio-quiet 
quasars, the open diamonds represent radio-loud quasars, and the asterisks represent NLSls. 
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Fig. 2. - Composite spectra, sorted by L/L^^^, with the mean L/L^^^ values indicated at the right. The top spectrum is a composite 
spectrum created only from objects identified in the literature as NLSls. The horizontal dashed lines and tick marks indicate the normalized 
continuum levels for each composite spectrum. The scaled height of the normalized continua above zero flux, as indicated for the bottom 
spectrum, is the same for all spectra plotted. 
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Fig. 3. — Multi-component Gaussian fits (smooth solid curves) to observed emission lines (jagged solid curves) for a composite spectrum 
of < Zz/I/edd >~ 0.58. The continuum is normalized by a power law continuum fit. The continuum (at unity) and composite line fits are 
plotted. The individual components of the fits are shown for the fit to N III] as dotted lines. In the other panels, the dotted lines represent 
the sum contribution of each labelled emission line. There is an unidentified bump at ~1070 A that is not attributed to O VI (see Hamann 
et al. 1998). 
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Fig. 4. — Distribution of estimated Eddington ratio, L/Lg^d- The two arrows are the median (0.59) and mean (0.89) i/Ledd of the entire 
sample. About 27% of the objects have L/Ledd > 1- 
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Fig. 8. — Normalized composite spectra sorted by L/Lgdd (top left, identical to Fig. 2 without the NLSl composite) axe compared to 
normalized composite spectra computed from different ranges in bolometric luminosity (top right, in units of ergs/sec), SMBH mass (bottom 
left, in M0), and FWHM(C IV) (bottom right, in km/sec). Dashed hues and tick marks as in Fig. 2. The L, SMBH mass, and FWHM 
composite spectra all clearly show trends for decreasing line equivalent widths and peak heights as luminosity, SMBH mass, and FWHM(C IV) 
increase. The composite spectra sorted by FWHM(C IV) even exhibit this trend despite spanning less than half an order of magnitude in 
average luminosity, a much narrower range than in the L/Lgdd composites (see Table 1). However, the L/L^^^ composites show nearly 
constant peak heights and decreasing FWHMs with increasing L/Lg^^. 
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Fig. 9. — Composite spectra sorted by L/Ledd for a narrow range in SMBH mass of 10**Mq < Msmbh < IO^Mq (left) are compared to 
ones sorted by L/L^^id for narrow range in luminosity of 10^^ ergs/s < L < 10^** ergs/s (right). The spectra are shown prior to normalization 
to the eontiinium. The eoniposite spectra at nearly constant A/smbh show the same emission-line behavior as the L/Lgdd composites created 
from our entire sample: nearly constant peak heights and decreasing FWHMs with increasing L/L^^^^. In contrast, the composite spectra at 
nearly constant L do not share this behavior and show a trend similar to the Baldwin Effect: decreasing peaJj heights and equivalent widths 
with increasing Msmbh- 
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Fig. 10. — Metallicities derived from comparisons of diflferent line ratio diagnostics to theoretical results from Figure 5 in Hamann et al. 
(2002) are shown as a function of L/L^m- The average includes N III]/C III], N III]/0 III], and N V/(C IV+O VI). Symbols as in Fig. 7. 
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Fig. 11. — Symbols as in Fig. 7. Our best estimates of the overall metallicity of each spectrum are shown for the L/L^^^, SMBH mass, 
and luminosity composite spectra. The NLSl composite, plotted in all three panels as a filled triangle, has a metallicity that is slightly high 
for its SMBH mass and luminosity, but still well below the metallicities of most luminous quasars. 



Table 1 
Composite Parameters 



L/Ledd a # Objects z Log L FWHM(C IV) Log Msmbh 
at C IV [ergs s'^] [km s^^] [Mq] 



L/Ledd Composites 



0.17 


-0.65 


64 


1.06 


46.84 


6900 


9.48 


0.30 


-0.56 


62 


1.17 


47.17 


5600 


9.60 


0.41 


-0.49 


106 


1.62 


47.16 


5100 


9.44 


0.58 


-0.40 


85 


1.83 


47.43 


4600 


9.56 


0.83 


-0.60 


90 


1.93 


47.37 


3900 


9.36 


1.38 


-0.33 


110 


2.55 


47.63 


3400 


9.42 


3.41 


-0.54 


45 


2.86 


47.84 


2400 


9.26 



NLSls 

0.67 -0.93 18 0.06 45.53 2900 7.83 



L Composites 



0.43 


-0.73 


15 


0.57 


44.73 


3100 


7.25 


0.46 


-0.66 


90 


0.37 


45.69 


3800 


8.10 


0.69 


-0.36 


155 


1.39 


46.62 


4400 


8.87 


1.15 


-0.47 


278 


2.63 


47.52 


4900 


9.59 


1.21 


-0.24 


20 


3.01 


48.32 


5900 


10.30 



-^SMBH Composites 



0.61 


-0.98 


7 


0.56 


44.62 


2500 


6.78 


0.69 


-0.90 


61 


0.46 


45.68 


3000 


7.27 


1.08 


-0.71 


198 


1.43 


46.84 


3800 


8.65 


0.84 


-0.59 


261 


2.46 


47.52 


5200 


9.54 


0.58 


-0.57 


34 


2.84 


48.17 


7500 


10.28 








FWHM(C IV) Composites 
















3.94 


-0.91 


14 


1.64 


47.25 


1700 


8.35 


1.24 


-0.48 


237 


1.71 


47.22 


3100 


8.87 


0.58 


-0.38 


195 


1.86 


47.48 


4900 


9.43 


0.38 


-0.42 


91 


2.29 


47.65 


6700 


9.84 


0.17 


-0.47 


25 


1.89 


47.52 


9300 


10.07 



Table 2 



Emission Line Data 





Property 


O VI A1034 


Lya A1216 


N V A1240 


C IV A1549 


He II A1640 


O III] A1665 


N III] A1750 


C III] A1909 


0.17 


Flux/Lya 


0.21 


1.00 


0.32 


0.75 


0.10 


0.07 


0.06 


0.22 




REW" 


16 


80 


26 


80 


12 


9 


7 


32 




FWHM'' 


5500 


4700 


6900 


6800 


6800 


6800 


6900 


4900 


0.30 


Flux/Lya 


0.17 


1.00 


0.21 


0.57 


0.08 


0.06 


0.04 


0.13 




REW 


16 


94 


20 


74 


12 


8 


6 


23 




FWEM*" 


5200 


4300 


5100 


4800 


4800 


4800 


5000 


5000 


0.41 


Flux/ Ly Of 


0.21 


1.00 


0.23 


0.57 


0.08 


0.06 


0.05 


0.11 




REW 


19 


95 


22 


76 


12 


9 


7 


20 




FWHM'' 


6200 


4600 


5300 


4900 


4900 


4900 


5100 


4300 


0.58 


Flux/Lya 


0.28 


1.00 


0.24 


0.57 


0.08 


0.06 


0.05 


0.14 




REW" 


22 


83 


20 


67 


10 


8 


7 


23 




FWHM'' 


6500 


3800 


4300 


4000 


4000 


4000 


4200 


4900 


0.83 


Flux/Lya 


0.22 


1.000 


0.19 


0.50 


0.07 


0.06 


0.03 


0.12 




REW" 


20 


96 


18 


65 


10 


9 


5 


20 




FWHM'' 


5200 


3000 


4300 


3700 


3700 


3700 


4000 


4000 


1.38 


Flux/Lya 


0.26 


1.00 


0.21 


0.47 


0.07 


0.05 


0.04 


0.11 




REW" 


21 


84 


18 


57 


10 


7 


6 


19 




FWHM* 


5500 


2800 


3900 


3300 


3300 


3300 


3600 


3500 


3.41 


Flux/Lya 


0.21 


1.00 


0.15 


0.41 


0.06 


0.04 


0.04 


0.10 




REW" 


20 


100 


15 


57 


9 


7 


7 


19 




FWHM* 


5100 


2100 


3000 


2500 


2500 


2500 


2800 


3100 


0.67 


Flux/Lya 




1.00 


0.14 


0.29 


0.06 


0.03 


0.03 


0.13 


(NLSls) 


REW" 




122 


18 


45 


10 


4 


5 


25 




FWllM'' 




2100 


281K) 


2200 


2200 


2200 


2700 


:!1()() 



"In units of A 
''In units of km s 



